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1. Introduction 


Abstract 

Objectives: The outbreak of Middle Eastern respiratory syndrome coronavirus 
(MERS-CoV) was one of the major events in South Korea in 2015. In particular, this 
study pays attention to formulating a mathematical model for MERS transmission 
dynamics and estimating transmission rates. 

Methods: Incidence data of MERS-CoV from the government authority was 
analyzed for the first aim and a mathematical model was built and analyzed for 
the second aim of the study. A mathematical model for MERS-CoV transmission 
dynamics is used to estimate the transmission rates in two periods due to the 
implementation of intensive interventions. 

Results: Using the estimates of the transmission rates, the basic reproduction 
number was estimated in two periods. Due to the superspreader, the basic 
reproduction number was very large in the first period; however, the basic 
reproduction number of the second period has reduced significantly after 
intensive interventions. 

Conclusion: It turned out to be the intensive isolation and quarantine in- 
terventions that were the most critical factors that prevented the spread of the 
MERS outbreak. The results are expected to be useful to devise more efficient 
intervention strategies in the future. 


MERS-CoV exhibited outbreaks in several regions and 
the last one was in South Korea (at the time of writing). 


One of the most significant events in South Korea in __ It is known as having a high fatality rate and having a 
2015 was the introduction and diffusion of Middle resemblance in clinical features with severe acute res- 
Eastern respiratory syndrome coronavirus (MERS- piratory syndrome coronavirus (SARS-CoV) [2], thus, it 
CoV). It is a novel coronavirus which was first identified has drawn much concern both from the public health 
in Saudi Arabia in September 2012 [1]. After that, authority and the general public when the first index 
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case was identified in May 2015, for the first time in 
Korea. 

In the case of the outbreak in South Korea, it is 
characteristic in that a single importation of the path- 
ogen sparked a large cluster of infected cases, while 
sporadic importations did not lead to local outbreaks in 
other regions [3]. It was surprising because South Korea 
has been said to have one of the most advanced medical 
and public health systems in the world [4], and it raises 
fears that the situation would be much worse in coun- 
tries with less developed public health systems. Thus, 
the investigation of this study about the South Korean 
incidence is expected to give an invaluable lesson both 
to public health authorities around the world and the 
community of epidemiology researchers. Moreover, it 
might be able to provide a quantified result with math- 
ematical modeling and statistical analysis, which is 
beneficial given that the studies about the South Korean 
case are largely comprised of epidemiological de- 
scriptions [4,5]. 

It is not only in the South Korean case but also in 
MERS-CoV studies in general that case studies and 
epidemiological descriptions are the main approach. Assiri 
et al [6] provided a detailed description of cases in the 
Middle East, Drosten et al [7] and Guery et al [8] analyzed 
the clinical features of infected cases, and Memish et al [9] 
and Al-Tawfigq et al [1] described the epidemiological data 
with respect to family clusters and hospitalized patient 
respectively. The most thorough description of the South 
Korean case is presented by Korea Centers for Disease 
Control and Prevention [5], while it was briefly mentioned 
by Chowell et al [3] and Cho and Chu [4]. 

Other than these descriptive studies, a small number 
of literature tried to assess the transmissibility of MERS- 
CoV in a more quantified fashion. Cachemez et al [10] 
estimated the incubation period and generation time 
from publicly available data and calculated the repro- 
duction numbers both in animal-to-human and human- 
to-human transmission cases. Breban et al [11] 
compared the possibility of MERS-CoV being 
pandemic by comparing the basic reproduction number 
(Appendix) with that of SARS. Using MERS-CoV 
outbreak data, Chowell et al [12] delved further and 
calculated two reproduction numbers based on different 
epidemiological scenarios: one is when all reported 
zoonotic cases are severe and the other is not. Chowell 
et al [3] took a different approach that compared the 
reproduction numbers in SARS and MERS-CoV in large 
hospital clusters. 

Since South Korea is where the largest number of cases 
was reported outside the Middle East, we need to pay more 
attention to the South Korean case. While a few researches 
have dealt with the South Korean case, as briefly 
mentioned above, research gaps still exist that need to be 
filled. First of all, whether some demographic character- 
istics, such as age or sex, and locations have influence on 
the possibility of being infected needs to be explored. Also, 


Y. Kim, et al 


the transmissibility needs to be investigated with more 
detail. As done in some literature, transmissibility in 
hospitals and in general need to be calculated separately 
and it needs to be taken into consideration whether the 
transmissibility might change in time. Here, the SIR-type 
mathematical modeling is of great use, while there are very 
few in MERS-CoV studies. Its usefulness might be boos- 
ted if combined with real-world data, which will be done in 
the present study. 

The purpose of this study is to uncover the charac- 
teristics of MERS-CoV transmission in South Korea in 
2015. To this aim, a mathematical model will be built 
and analyzed. In particular, it will be observed whether 
there were any periods with different transmissibility 
and, if so, they will be estimated in general and hospital 
circumstances. Before this, some statistical analysis 
will be conducted in order to identify which group was 
more susceptible to MERS-CoV based on South 
Korean data. 


2. Materials and methods 


A mathematical model of MERS-CoV transmission 
in South Korea is suggested based on the model in 
Chowell et al [12]. It categorizes each individual into 
one of six epidemiological classes; susceptible (S ), 
exposed (or high-risk latent) (E ), symptomatic and in- 
fectious (J), infection but asymptomatic class (A), 
hospitalized (H ), and recovery class (R). It is assumed 
that only infectious and hospitalized individuals can 
infect others and asymptomatic individuals cannot. In 
Chowell et al [12], the actual data of the zoonotic case 
were gathered so they were able to take secondary cases 
as well as index cases into account. By contrast, there 
was no zoonotic transmission except the primary case in 
South Korea, thus we do not consider the zoonotic case 
of MERS-CoV model in the present study. The model 
takes the following form: 


dS _ s BU +/1H) 

dt N 

dE _ 6(1+1H) 

7 S N KE 

dl 

a KOE > (Ya ar 1) I 
(1) 

dA 

—=k(l1—p)E 

HK 8) 

dH 

—_—=y,l-y,H 

at Ya Yr 

oi I++7,H 

dt =71 Yr 


where 6 is the human-to-human transmission rate per 
unit time (day) and / quantifies the relative trans- 
missibility of hospitalized patients; x is the rate at which 
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an individual leaves the exposed class by becoming 
infectious (symptomatic or asymptomatic); p is the 
proportion of progression from exposed class E to 
symptomatic infectious class J, and (1 — p) is that of 
progression to asymptomatic class; A; y, is the average 
rate at which symptomatic individuals hospitalize and y, 
is the recovery rate without being hospitalized; y,. isthe 
recovery rate of hospitalized patients. Using these 
parameters, we will estimate the transmission rate (() 
and transmissibility of hospitalized (/ ) in next section. 

In many epidemiological models, basic reproduction 
number is one of the key values that can predict whether 
the infectious disease will spread into a population or 
die out. It is the average number of secondary infectious 
cases when one infectious individual is introduced in a 
whole susceptible population. The basic reproduction 
number is calculated by the next generation matrix 
approach outlined in van den Driessche and Watmough 
[13]. As a result, the basic reproductive number Ro for 
our model is: 


Ro= pb ( 1, __Yel ) (2) 


Vor, What yD 


where the first term is the average number of secondary 
infected people from one infectious individual during 
his or her infectious period and the second term is the 
average number of secondary infected people from one 
patient during their hospital stay (the details of the 
derivation of Ro are given in Appendix). 


3. Results 


3.1. Identifying more susceptible groups 

Other than calculating transmissibility, revealing 
more susceptible groups is one of the useful ways of 
characterizing the transmission of MERS-CoV in South 
Korea. It is done by conducting statistical analysis to 
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publicly available MERS-CoV data in South Korea and 
its result is presented in Table | (this is based on Tables 
3 and 4 and Figures 2—4). 

First of all, it was tested whether age is related to the 
susceptibility and mortality of MERS-CoV using z-test. 
The average age of MERS-CoV confirmed cases is 
55 years and that of death patients is 69.8 years, while 
that of the total population in South Korea is 38.9 years 
in 2015. These differences were statistically significant 
and age is related both to susceptibility and mortality of 
MERS-CoV. Next, Chi-square test was conducted in 
order to investigate whether sex, region, and/or visited 
hospital had any relationship with mortality of MERS- 
CoV. The result shows that mortality was different by 
region and visited hospital, while it was not by sex. We 
also tested whether patients’ sex was different by age, 
region, and/or category of patients (health workers, 
epidemiological survey, or visitor/family). According to 
the result, the sex of patients was not related to age and 
region, but it was significantly different by the category 
of patients. Finally, we investigated whether sex, region, 
and mortality are related to the number of cases by 
generation; in other words, it was investigated whether 
mortality increased or decreased as the disease spread 
further. The result of Chi-square test demonstrates that 
the number of cases by generation was not related to age 
while it was different by region and mortality (more 
details are given in Appendix). 


3.2. Estimating transmission rates 

The parameter estimation of system (1) to the inci- 
dence data of MERS-CoV in South Korea from May 6, 
2015 to July 1, 2015 is conducted using the least-square 
method for this nonlinear system. Since the primary case 
was exposed to MERS-CoV during their stay in the 
Middle East, the patient and his wife are assumed to be 
in the exposed class of our system before the confir- 
mation day (May 20, 2015). 


Table 1. Results of statistical analysis of Middle Eastern respiratory syndrome coronavirus data. 


Susceptibility 
IV Age z= 12.9 
(p < 0.001)* 
Sex — 
Region = 
Mortality — 
Category = 


Visited hospital 


*Significant at 0.05 level; ‘Significant at 0.1 level. 


DV 
No. of cases 
Mortality Sex by generation 
Z = 8.65 x = 62.1158 - 
(p < 0.001)* (p = 0.472) 
x2 = 0.5818 — ¥° = 0.9321 
(p = 0.4456) (p = 0.6275) 
° = 11.1589 7° = 6.3462 ° = 52.9542 
(p < 0.001)* (p = 0.1747) (p < 0.001)* 
— _ = 4.591 
: (p = 0.1007)' 
— x? = 12.6924 — 
(p < 0.001)* 
y° = 13.5511 — _ 
(p < 0.001)* 
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Figure 1. Data and parameter estimation of two periods. 


In Figure 1, the data of MERS-CoV incidence (°c) 
and its best fit (solid curve) of incidence Eq. (1) are 
plotted from May 6, 2015 to July 4, 2015. The MERS- 
CoV incidence is defined as the number of newly 
infected individuals per day. The figure demonstrates 
that our fitted MERS-CoV model fits well to the re- 
ported incidence data. In the figure, we can see the trend 
of logistic curve; thus, we split up the time into two 
periods (May 6, June 7 and June 7, July 1) where the rate 
of change is in its increasing and decreasing trend 
respectively. Due to the implementation of intensive 
interventions at June 7, each period has its own 
remarkable trend, thus, parameters are estimated in each 
period. Two kinds of transmission rates are estima- 
ted—transmission rate 6 and transmissibility of hospi- 
talized patients /—because many MERS-CoV patients 
are infected in hospitals and it is assumed that the value 
of hospitalized transmissibility / is bigger than trans- 
mission rate 6. Parameter values used for simulation are 
listed in Table 2. 

As a result, the estimated value of the transmission 
rate 6 is 0.0835 [95% confidence interval (CI): 
0.0835—0.0835] and the transmissibility of hospitalized 
patients / is 22 (95% CI, 18.7303—25.2697) in Period 1, 


Table 2. Parameter used in simulations for Figure 1. 


Parameter Description Value Refs 

B human-to-human Estimated 
transmission rate 

I transmissibility of Estimated 
hospitalized patients 

K leaving rate from E class 1/(6.6) [14] 

p fraction from E' to J 0.585 [13] 

We average rate from / to H 0.6403 [15] 

Yr average recover rate without 1/5 [15] 
being hospitalized in J 

Yr Rate from H to R 1/7 [15] 
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and 6 is 0.036 (95% CI, 0.0306—0.0414) and /is 1 (95% 
CI, 0.1597—1.8403) in Period 2. Also, as mentioned in 
the previous section, the basic reproduction number is 
calculated. Using our estimated values of 6 and /, Ro is 
5.3973(95% CI, 4.6057—6.1888) in Period 1 and 0.1351 
(95% CI, 0.0403—0.2562) in Period 2. The remarkably 
larger value of Ro in Period 1 than in Period 2 is 
considered to be attributed to the superspreader. 


4. Discussion 


Due to the rapid development of nation-wide and 
world-wide transportation, it is more likely for a large- 
scale infectious disease outbreak to occur. This raises 
the importance of timely and accurate countermeasures. 
In order to devise efficient intervention strategies, it is 
critical to understand both the biological characteristics 
of the disease and the social response of the location 
where the disease spreads or is likely to spread. It is 
because these two aspects in combination determine the 
mode of diffusion of the infectious disease. In this re- 
gard, it is crucial to investigate the characteristics of 
MERS-CoV transmission in South Korea since it can 
demonstrate the current capability and the future po- 
tential of South Korea to resist infectious disease. This 
study tried to meet this need by examining which de- 
mographic group is more susceptible using incidence 
data and estimating transmission rates of general and 
hospitalized populations. The result of this study is ex- 
pected to have theoretical and practical implications. 

Further research can be suggested as follows. As seen 
in the result, the role of the superspreader might be 
decisive in the spread of MERS-CoV in South Korea. It 
can be explored what might happen if the superspreader 
was identified and isolated earlier or later, or if the 
spreading power was stronger or weaker. Also, the role 
of the media in the spread of MERS-CoV in South 
Korea can be investigated. How the Korean media 
covered MERS-CoV and how the lay person responded 
can be an issue, and how these two combined influence 
the spread of MERS-CoV may be examined. 
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Appendix. The basic reproduction number 


The basic reproductive number is calculated by using 
the methodology (the next generation matrix approach). 
Now, we let F(x) represent the rate of appearance of 
new infections. The net transition rates out of the cor- 
responding compartment are represented by V(x). Then, 
we find the Jacobian matrix of F(x), V(x), and denote 
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F = [0F/0x;] $ and V =[0V/0x;] evaluated at disease 
free equilibrium point x*, which consists of S = N and 
the rest of them is zero. 
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returning on May 4, 2015 from the Middle East (April 
24""—May 4"). Further epidemiological investigation 
showed that the primary case had also traveled to the 
Barain, United Arab Emirates, Saudi Arabia, and Qatar. 


oo : 4 . oo - id : The first symptom, which included fever, muscle aches, 

00 0 ; . a M cough, and dyspnea, was detected on May 11", 2015. 

ae On May 20", 2015, he was diagnosed with MERS-CoV, 
pp x PY abl 6B abl pl 
tet Yat) Yat tr VYat 1) 
= 0 0 0 
0 0 0 


The matrix FV~! has two zero eigenvalues and one 
positive eigenvalue. 


1 Val 
Ro = vo | . ) 
. VatY1 Ve(Yat¥7) 


Description of MERS-CoV data 


We retrieved data for MERS-CoV in the Republic of 
Korea from public reported sources (the Korea Centers 
for Diseases Control and Prevention and the Korean 
Ministry of Health and Welfare). As of July 7, 2015, 186 
laboratory-confirmed MERS-CoV cases and 36 fatal 
cases have been reported. 

The earliest two confirmed cases were the primary 
case (68-year-old man) and his wife. The ongoing 
outbreak in South Korea began when the primary case 
developed respiratory illness on May 11, 2015 after 


Table Al. Characteristics of confirmed cases. 

All cases Fatal cases 
Characteristics (n = 186) (n = 36) 
Age group (y) 
0-18 1 0 
19-39 39 0 
40—59 71 1 
60—79 66 22 
> 80 9 7 
Sex 
Male 111 24 
Female is 12 
Occupation 
Health workers 39 2) 
Epidemiological survey 82 23 
Visitor/family 65 11 


so the patient was isolated and treated in hospital 
designated by the Korea government. After that, MERS- 
CoV spread in South Korea from the man, and we 
currently have 186 confirmed cases including 36 death 
cases (July 7, 2015). 

Table Al shows three kinds of characteristics of 
confirmed cases; age group, sex, and occupation. In 
occupation, infected patients are classified into three 
groups; Health Workers, Epidemiologic Survey, and 
Visitor/Family. We can derive the mortality rates by 
occupation from the numbers in the table, which are 
2.6% in health workers, 25.6% in Epidemiologic Sur- 
vey, and 17% in Visitor/Family, and we can see that the 
mortality rate is the highest in Epidemiologic Survey. 

Figure Al presents the time, the location, and the 
number of infections by the primary patient until he was 
isolated, and Table A2 displays the number of infected 
cases by the hospitals that the primary patient had 
visited before isolation. As the table shows, while it was 
Pyeongtaek St. Mary’s Hospital where the primary case 
was confirmed, Samsung Medical Center was the loca- 
tion where the largest number of people were infected. It 
also draws attention that the number of cases decreases 
by the order of infection in Pyeongtaek St. Mary’s 
Hospital while it tends to increase in Samsung Medical 
Center. 

Figure A2 illustrates a transmission tree that MERS- 
CoV spread from primary case to 153 cases on June 15, 
2015. The number of first-order infected cases is 39, that 
of second-order infected cases is 104, and that of third- 
order infected cases is seven. We omit some patients in 
Figure A2 since the infection route of Patient No. 119 is 
unknown. Patient No. 126 stayed in the same location 
with Patients No. 14, No. 17, No. 26, and No. 108, so the 
infection of this patient cannot be attributed to one other 
infector. As the figure indicates, Patients No. 1, No. 14, 
and No. 16 have infected considerably more patients 
(39, 76, and 21 people, respectively) than usual cases, 
and we can see clearly that the MERS-CoV spread in 
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Primary case 


Pyeongtaek St. Mary's H 
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The first symptoms 
occurrence 


* May 12, May 14, 
May 15 (outpatient) 


— 1 people infection 


Figure Al. 


Table A2. The number of infections by the hospital that 
the primary patient visited. 


Pyeongtaek Samsung 
Asan Seoul 365 Medical St. Mary’s Medical 


Clinic Center Hospital Center 
Secondary 1 1 26 0 
a 1 0 10 74 
4 0 0 0 3 
South Korea is largely attributed to these 


superspreaders. 

The number of patients of MERS-Cov on a daily 
basis from May 20, 2015 to July 1, 2015 and the hospital 
where those patients were confirmed are presented in 


Primary case _ 


Figure A2. 


Samsung Medical Center 


UY 


* Pyeongtaek-sungmo Hospital 
: May 15 — May 17 (hospitalization) 


* 365 medical : May 17 (outpatient) 


* May 17 (outpatient) 
* May 18 — May 20 (hospitalization) 


— 36 people infection — 89 people infection 


* May 18 Emergency observation 


= 1 people infection room isolation. 


The timeline of the primary patient before isolation. 


Figure A3. As the figure illustrates, the number of pa- 
tients reached its peak by 23 on June 7th, 2015. The 
figure also illustrates that the hospital where the largest 
number of patients were confirmed changed from Hos- 
pital B (Pyeongtaek St. Mary’s Hospital) to Hospital D 
(Samsung Medical Center) by the peak day. 

When it comes to the region where confirmed cases 
were infected, Gyeonggi-do, Seoul, and Daejeon are the 
three major infection areas where their total number of 
infected people is 147 which accounts for about 80% of 
the total number of infected people, and Seoul metro- 
politan area is the region where the majority of in- 
fections occurred (figure omitted). In each region, the 
hospitals where the most infections occurred are as 
follows: Samsung Medical Center in Seoul, Pyeongtaek 
St. Mary’s Hospital in Gyeonggi-do, and Dae-Cheong 


Transmission tree of confirmed Middle Eastern respiratory syndrome coronavirus. 


Figure A3. 
syndrome coronavirus. 


Middle Eastern respiratory syndrome coronavirus 
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Incidence data of MERS-CoV, South Korea 
May 20, 2015 to July 1, 2015 (n = 186) 


N 
ul 


01/Jun 
28 

24 
20/May 


® Hospital A 
B Hospital G 


B Hospital B 
®@ Hospital H 


Hospital and Konyang University Hospital in Daejeon. 
The numbers of hospital infections are 89, 36, 14, and 
11, respectively. The combined number of infection in 
other hospitals is < 10. 
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